On the saturation of the centrifugally excited 
curvature drift instability in AGN magnetospheres 



Osmanov Z. 

Georgian National Astrophysical Observatory, Ilia Chavchavadze State 
University 

Kazbegi ave. 2", Tbilisi, 0160, Georgia 
E-mail: z.osmanov@astro-ge.org 

Abstract. We study a saturation process of the centrifugally driven curvature 
drift instability (CDI) in AGN magnetospheres close to the light cylinder surface 
to examine the twisting of magnetic field lines leading to the free motion of AGN 
winds, that completely kills the instability. Considering the Euler, continuity, and 
induction equations, by taking into account the resonant conditions, we derive 
the growth rate of the CDI. We show that due to the centrifugal effects, the 
rotational energy is efficiently pumped directly into the drift modes, that leads 
to the generation of a toroidal component of the magnetic field. As a result, the 
magnetic field lines transform into such a configuration when particles do not 
experience any forces and since the instability is centrifugally driven, at this stage 
the CDI is suspended. 
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1. Introduction 

One of the fundamental problems concerning AGN relates to the understanding of a 
question: how the plasma goes through the light cylinder surface LCS (a hypothetical 
zone where the linear velocity of rigid rotation exactly equals the speed of light). It is 
generally believed that in AGN magnetospheres the magnetic field is strong enough to 
provide the frozen-in condition of plasmas, flowing in the direction of corotating field 
lines. On the other hand, such a dynamics of plasma particles will inevitably lead 
to a situation when particles' velocity achieves the speed of light on the LCS, thus 
violating the causality principle. Therefore, it is clear that in the mentioned zone, a 
certain twisting process of the magnetic field lines must exist. In [TJ Rogava et al. 
considered curved trajectories and generalized a work developed by Machabeli and 
Rogava [2] ■ The authors examined the dynamics of a single particle, that moves along 
a corotating, curved channel. In this simple mechanical model the channel plays a role 
of magnetic field lines. It was shown that the dynamics of particles asymptotically 
becomes force-free if a shape of the channel is given by the Archimedes' spiral. It 
has been found that in the aforementioned case the particles cross the LCS without 
violating the causality principle. 

Generally speaking, if the magnetic field is still robust, one has to twist the 
field lines in an appropriate way. Therefore, it is necessary to generate the toroidal 
component of the magnetic field provided by certain current. It is well known that 
particles moving along curved field lines, also drift perpendicularly to the curvature 
plane even if the curvature is very small. Such a drift motion creates curvature current, 
which leads to the necessary condi- tions for the development of the CDI (see e.g., 

The innermost region of AGN magnetospheres is supposed to be rotating, 
therefore, the centrifugal force (CF) seems to be essential for studying the 
magnetospheric plasma motion. Magneto-centrifugal effects have been extensively 
studied in a series of papers. In [5] the angular momentum and energy pumping process 
from accretion disks have been considered. The authors emphasized the special role 
of CF in dynamical processes governing the acceleration of plasmas. The role of the 
centrifugal acceleration in producing the non-thermal radiation from rotating AGN 
winds has been studied in [6]. We generalized this work in [7] and have shown that 
under certain conditions electrons might reach very high Lorentz factors 10 5-8 . The 
similar investigation performed for non-blazar objects was developed in [8] and the 
high efficiency of centrifugal acceleration was emphasized. 

Generally speaking, the centrifugal acceleration may induce plasma instabilities. 
When the CF acts on a particle and changes with time, it plays a role of a parameter, 
and drives the so-called parametric instability. Ccntrifugally excited unstable plasma 
waves deserve a special interest in different astrophysical scenarios. In [9] we considered 
the corotating Crab magnetspheres arguing that the centrifugal force inevitably causes 
the separation of charges, that in turn leads to the creation of the Langmuir waves. 
We have shown that due to the centrifugal effects the instability is very efficient. This 
method was applied to AGN jets in [10] where we studied the stability problem of 
the rotation-induced electrostatic waves. The CF is significant also also for inducing 
the CDI. Even if the field lines have a very small curvature, it might cause a 
drifting process of plasma, creating current, which will inevitably produce the toroidal 
component of the magnetic field. However, a value of the mentioned component is 
not strong enough to change a field line's configuration significantly. In this context 
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corotation plays an important role. Namely, as is shown in our model by means 
of the centrifugal acceleration, the rotational energy can be pumped directly into the 
toroidal component, amplifying it efficiently. This is the meaning of the curvature drift 
instability. For examining the role of the corotation in the curvature drift instability 
for pulsar magnetospheres, in [11] we studied the two-component relativistic plasma, 
estimating the increment of the instability. We have found that the growth rate was 
bigger than pulsar's spin-down rates by many orders of magnitude, indicating high 
efficiency of the CDI. This instability is very important since, due to the mentioned 
drift, current is produced, which, via the parametric mechanism leads to the creation 
of the amplifying toroidal magnetic field, reconstructing the magnetosphere. By the 
influence of the toroidal magnetic field, the initial field lines are twisting until they 
transform into the shape of the Archimedes spiral, when the motion of the particles 
is described by the so-called force- free regime [12] , 

In the present paper we investigate the saturation (transition) mechanism of the 
CDI in AGN winds. 

The paper is arranged as follows. In Sect. 2, we introduce the curvature 
drift waves and derive the dispersion relation and a corresponding expression of the 
transition timescale. In Sect. 3, the results for typical AGN are presented and, in 
Sect. 4, we summarize our results. 



2. Excitation of curvature drift waves 



We consider plasma that consists of relativistic electrons with the Lorentz factor 
7 ~ 10 5-8 [7j, [5]. It is assumed that the field lines initially are almost rectilinear 
and corotate. As we have already mentioned in the introduction, current providing 
the twisting of field lines is created due to the curvature drift, characterized by the 
following (curvature drift) velocity 

lav 2 

"=T^> (1) 

where = eBg/mc, e and m are particle's charge and the rest mass respectively, Bq 
is the unperturbed magnetic induction, c is the speed of light, Rb is the curvature 
radius of magnetic field lines, 70 is the initial Lorentz factor of particles and is the 
longitudinal velocity. As is clear from Eq. Jl}, the drift velocity is proportional to 
7om, therefore, the corresponding value of the bulk flow (protons) (having the Lorentz 
factor of the order of 10), will be by many orders of magnitude less than that of the 
relativistic electrons (with 70 ~ 10 5-8 ). This in turn means that the contribution of 
protons is negligible and we consider one-component plasma composed of relativistic 
electrons. 

In the zeroth approximation particles move along the magnetic field lines and 
our aim is to consider the twisting process of these field lines and study a subsequent 
saturation mechanism. 

For studying the development of the CDI we consider the Euler equation, which 
governs the dynamics of corotating plasma particles [13j 

^ + (v • V)p = -c 2 7£V£ + - f E + -v x bV (2) 
at m \ c J 

where 



£, = v/l-fi 2 i? 2 /c 2 , 
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Figure 1. a) The geometry in which the set of main Eqs. i 12!4[ i is 
considered; (eg,e r ,e x ) denotes orthonormal basis of unit vectors; is directed 
perpendicularly to the plane of the figure (equatorial plane); C is the center of 
the curvature; and eg and e r are the tangential and perpendicular (with respect 
to the field line) unit vectors respectively, b) The geometry for deriving Eq. J26}: 
the curved line denotes the twisted magnetic field, B, generated due to the raising 
of magnetic perturbation, B r . Note, that B r and Bg are oriented with respect 
to the initial 'quasi straight' magnetic field line. 



and p is the momentum, v is the velocity, 7 is the Lorentz factor of the relativistic 
particles, and E and B are the electric field and the magnetic induction, respectively. 
fl is the angular velocity of rotation and R is the coordinate along the magnetic field 
lines. We express the equation of motion in the cylindrical coordinates (see Fig. [T^,). 
The first term on the right-hand side of the Euler equation — c 2 7£V£ represents the 
centrifugal force. As we see, this force becomes infinity on the LCS, therefore, its 
overall effect is significant in the mentioned zone. For closing the system, we add to 
Eq. ^ the continuity equation: 

| + V.J = 0, (3) 

and the induction equation: 

„ 13E 4i 

VxB = -— + — J, 4 

c at c 

respectively. By p = en and J = env we denote the charge density and the current 
density, respectively and n represents the electron density. 

If we take into account the frozen-in condition, E + v x B /c = 0, describing 
the leading state of the system, then one can show that Eq. @ reduces to [5] : 



dv n 2 R 



dt 1-0^ 



fl 2 R 2 2v 2 



(•5) 



For the ultra relativistic case (7 >> 1) and the following initial conditions: R(Q) = 0, 
v(0) w c (v = dR/dt) a solution to Eq. (O is given by 

v(t) = v.. » ccos(fii + ip), (6) 

where (p denotes the initial phase of a particle. 

As we have already mentioned above, we assume that the magnetic field lines 
initially have a small curvature, that drives the particles along the x axis (see Fig. 
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The mentioned drift causes current, which in turn creates the toroidal component 
of the magnetic field. 

For studying the development of the CDI we linearize the system of equations in 
Eqs. J2H1]), by expanding all physical quantities around the leading state 

+ (7) 

* = {n,v,p,E,B}, (8) 

where the zeroth order quantities are given by no = const, vq x = u, vqq — o. , 
Po = v /7o, B oe = B (see below equation (US])), B 0r = B 0x = 0, E = -v x B /c 
We then express the perturbed quantities as 

^\t,r) cx * x (t)exp[i(kr)] , (9) 

which reduces the set of equations in Eqs. (J2EJ) to 

% -i(k x u + k e v „)p£ = — v^Bl, (10) 
at " mc " 

- i{k x u + kev l{ )n 1 = ik x n°vl, (11) 

- ik e cBl = 4ire(n vl + n 1 ^), (12) 

where k x and kg are the wave vector components corresponding to the coordinate 
system described in Fig. ((T^,). By 



2L 

we denote the equipartition magnetic induction on the LCS for the leading state, 
(where L is the luminosity of AGN and R c = c/Q is the light cylinder radius). In 
deriving Eqs. (|10II12[) , the wave propagating almost perpendicular to the equatorial 
plane, was considered and the expression ~ cEl/Bg was taken into account. For 
simplicity, the set of equations are given in terms of the coordinates of the field line 
(see Fig. EJi). 

Introducing a special ansatz for v x and n 1 

V i = I4e lkA(t) , (14) 

n^ne^W, (15) 

where 

ut u 

A ^ = Y + 4n s[li[2{nt + v)] ' (16) 

A+(t) = ^sm(nt + cp), (17) 

and by substituting Eqs. (fl"4|) and (fl"5|) into Eqs. (|10H12|1 . it is straightforward to solve 
the system for the toroidal component 

-ik^cBUt) = ^e ikA W f e- lkA ^vJt')B r (t')dt'+ 

7e0 J 



k x ue tkA ^ J dt' J e-* A{t " ) v n (t")B r (t")dt" 
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(18) 

where w e = e\J Aim^/m and n e Q are the plasma frequency and initial electron density 
respectively. To simplify Eq. ([TS]) , one may use the following identity 



^iiix sin y 



^2 j s ( x )< 



±isy 



(19) 



where J s (x) (s = 0; ±1; ±2 . . .) is the Bessel function of integer order [14]. Eq. ([18 
then reduces to 



B r (uo) 



J s (g)J q (h)Ji(g)J P (h) 
2-feokec ^ ^ to + %H + ft(2s + q) 

a— ±1 s,q,l,p 2, V ^/ 



E E 



xB r (w + Jl(2[s-l] +<?-p + cr)) 



1 - 



+ ft (2s + q) 



QJ 2 e k x u 



E E 



Xe «?(2[s-(]+g-p+a-) I 

J.(g)J q (h)Ji(g)J p {h) 



-x 



' ^=±i », 9 ,i,p ( w + V + fi ( 2 [ s + A*] + 
x S r (w + ft (2[s - Z + /x] + g - p + a)) x e ^( 2 [ s -'+^+9-P+ CT ) , (20) 



where 



k x u k e c 

5 = lft' 

For solving Eq. (|20[) we must examine similar expressions, by rewriting Eq. (|20[) 
for the infinite number of components B r (uj ± ft), B r (u> ± 20),... etc. Therefore, 
the system becomes composed of the infinite number of equations, that makes the 
task unsolvable. To overcome this problem we use a physically reasonable cutoff on 
the-right hand side of the equation d5J. As is clear from Eq. ([20]) . the considered 
instability is characterized by proper frequency of the curvature drift modes 

w «-^p (21) 

when the corresponding conditions k x u/2 < 0, 2s + q = 0, and 2[s + fi] + q = are 
satisfied. As we see, only the resonance terms provide a significant contribution to the 
result. 

We consider the typical parameters for AGN magnetospheres L ~ 10 44 erg/s, 
ft = 3 x 10 _5 s _1 , 7 e o ~ 10 5 , Rb ~ -R; c n e o ~ 0.001cm -3 . Then, one can show 
that for the curvature drift waves with the wavelength A ~ Ri c (A = 2n/k. For 
making the physics realistic we have assumed that k x < and u > 0, otherwise the 
frequency becomes negative) one has the following relation \k x u/2\ ~ 10~ 12 s _1 << ft. 
Therefore, all terms with non-zero values of ft(2s + q) and ft(2[s + /i] + q) oscillate 
rapidly and therefore do not contribute to a final result. 

Since particles have different phases, to solve Eq. (|20|) . we must examine the 
average value of B r with respect to ip. Then, by taking into account the formula 

i- I e m *d^ = S N , , 
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Figure 2. Behaviour of /3r versus R/R c . The set of parameters is 0no 
{0.01; 0.3; 0.5; 0.7; 0.9}. 



after preserving only the leading terms the dispersion relation of the CDI is given by 

HQ 



ktrU 



uJ - 



Y-v J s(9)J q >(s^){h)Ji{g)Jp>(i,a)(h), 

cr./i— ±1 s.l 



where 



S = 2E±x = 



2^ eQ k e c 



(22) 
(23) 



a. 



To determine the CDI growth rate, T, let us write u> = uio + iT and substitute this 
into Eq. (f2"2"|) . Then, it is easy to show that the increment gets the following form 

2 

J2 J2 E ^9)J-2(s + ,)(h)Ji(g)J-2i + Ah) . (24) 

er,/x— ±1 s,l 

It is worth noting that this instability is unavoidable for plasma magnetospheric 
flows, because for the developing of the ICS a) the initial curvature should be 
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nonzero and b) the magnetic field must be robust enough to provide the frozen-in 
condition, (a)-is a necessary condition for creating the drift waves and (b)-guarantees 
the parametric mechanism of rotational energy pumping directly into the drift modes. 

3. Discussion 

Let us consider the Archimedes' spiral $ = aR, where $ and R are the polar 
coordinates and a = const. As was shown in !j) if the particle slides along the 
rotating channel that has the shape of the Archimedes' spiral then, an observer 
from the laboratory frame of reference will measure the effective angular velocity 
fl e f = fi + d$/dt = Q + av, where v is the longitudinal velocity of a particle. 
If the particle moves without acceleration in the laboratory frame of reference, the 
corresponding effective angular velocity must be equal to zero (£l e f — 0). Therefore, 
v — v c = —tt/a which means that, for any Archimedes' spiral with a < — f2/c, the 
dynamics of the particle becomes force-free and the particle moves with a certain 
"characteristic velocity", v c . 

After defining the corresponding components of the relativistic momentum Pr = 'ymv, 
Pq, = jmR£l e f, the equation of motion, dP/dt = F (F is the reaction force), leads to 
the following equation lj 



We note that (d 2 R/dt 2 ) = 0, when v = v c = — fi/a. 

Typical AGN outflows are highly relativistic, therefore, let us consider v c as c 
setting a = — fi/c. In Fig. [21 we plot the solution to Eq. ([25)) , that is the dependence 
of (3r = v/c versus R/R c (R c is the light cylinder radius) for different initial values of 
/3rq. As is clear from these plots, no matter what is the initial velocity of the particle, 
it asymptotically converges to the characteristic velocity, v c . This implies that the 
dynamics of the particle asymptotically becomes force-free. 

Now we can qualitatively analyze how the configuration of magnetic field lines 
changes with time. After perturbing the magnetic field in the transverse direction, the 
toroidal component will amplify. As a result the field lines will gradually lag behind 
the rotation. On the other hand, for the twisted field lines, their dynamical influence 
on particles will decrease. In the parallel regime, the efficiency of the CDI will also 
decrease and when the field lines get a shape of the Archimedes' spiral the instability 
completely vanishes since the CDI is centrifugally driven and as we see from Eq. (|25[) 
the radial acceleration becomes vanishing for this case. Therefore at this stage our 
aim is to estimate from Eq. ([24]) the timescale of transition of quasi-straight field lines 
into the Archimedes' spiral, in other words we have to estimate a characteristic time 
of saturation of the CDI. 

Let us assume that when the particle dynamics becomes force-free, the critical 
value of the toroidal magnetic field is B r , then referring to Fig. QJ>, one can show 
that tanx = B r /Bg. If we apply the exponential temporal behaviour of the toroidal 
component, B r ss i?J?exp(rr), then the corresponding transition timescale gets the 



d 2 R _ n-j 2 v{a + flv/c 2 ) 
dt 2 7 2 k 2 



ClefR, 



(25) 



where 
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Figure 3. The dependence of logarithm of the transition timescale on the 
normalized wavelength. The set of parameters is f! = 3 X 10 _5 s _1 , L/Le = 0.01, 
S e {10- 1 ;10- 2 ;10* 3 ;10- 4 }, 7 e0 = 10 s , R B ~ Rc, n e0 = 0.001cm" 3 and 
\ e = 2n/k = lOOfic 



(26) 



form 

r {b, 

We have taken into account that for the considered Archimedes' spiral (a = — O/c), 
tanx = 1. Generally speaking, the role of the CDI is twofold. At the one hand it 
guarantees the required twisting of magnetic field lines, and at the other hand, in 
terms of its feedback on plasma dynamics it provides necessary conditions for the 
saturation process. 

We consider the behaviour of the transition timescale versus the wavelength of 
the perturbation and the AGN bolometric luminosity respectively. For this purpose 
let us examine the typical AGN parameters: Mbh = 10 8 x Mq, fl = 3 x 10 _5 s^ 1 and 
L = I0 44 erg/s, where Mbh represents the AGN mass, Mq is the solar mass and L is 
the bolometric luminosity. 

In Fig. [3] we show the logarithm of the saturation timescale versus the wavelength 
normalized to the light cylinder radius for different values of initial perturbation 
B°/Bg = S e {1CT 1 ; 1CT 2 ; 1CT 3 ; 1(T 4 }. The set of parameters is fl = 3 x lO^s" 1 , 
L/L E = 0.01, 7 e o = 10 s , R B ~ R c , n e0 = 0.001cm- 3 and Xg = 2n/kg = 100i? c , where 
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Figure 4. The dependence of logarithm of the transition timescale on the 
dimensionless perturbation. The set of parameters is Q = 3 X 10~ 5 s~ 1 , L/Le = 



{0.01; 0.05; 0.3; 1}, y e0 
and X x = Rc 



10 s , R B « R c , n e0 = 0.001cm- J , \ g = 2-K/k e = 100R C 



Le = I0 46 erg/s is the Eddington luminosity for the given AGN mass. By different 
curves we show different cases of initial perturbation. As we see from the plots, the 
timescale is a continuously increasing function of X x (= 2ir/k x ). Indeed, from Eq. (|23|) 
we get 3o,±i ~ 1/A X , which combining with Eq. ([26]) provides r c< yA»- As a result, 
for the parameters considered in Fig. [31 the transition timescale varies from ~ 10 8 s 
{X x /R c = 1,8= lO^ 1 ) to ~ 10 9 s (X x /R c =2,6= 10- 4 ). 

In Fig. 31 we display the behaviour of log(r) versus the initial perturbation for 
different values of luminosities. The set of parameters is Q = 3 x 10 s , L/Le = 
{0.01; 0.05; 0.3; 1}, 7 e0 = 10 s , R B w R c , n e0 = 0.001cm- 3 , X g = 2n/k g = 100i? c and 
X x = R c - The figure shows the continuously decreasing behaviour of the transition 
timescale. This is a natural consequence of the fact that for bigger perturbations one 
needs lower time to reach the critical value B r . On the other hand as we see from 
the plots, the instability is less efficient (bigger timescale) for more luminous AGN. 
Such a behaviour is clearly seen from an expression of the drift velocity in Eq. |T]). 
Indeed, as we see, the drift velocity is proportional to the inverse value of the cyclotron 
frequency, which is bigger for bigger magnetic fields. On the other hand, from Eq. 
(|13p we have Bq oc \/L, which means that the curvature drift velocity is proportional 
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to i/vT. Since the curvature drift waves are more efficient for bigger curvature drift 
velocities, by increasing the bolometric luminosity, the increment of the instability will 
inevitably decrease and the corresponding transition timescale will be bigger. For the 
parameters used for Fig. Q] the transition timescale varies from ~ 10 10 s (S = 10 -4 , 
L/L E = 1) to ~ 10 8 s (S = lQ-\ L/L E = 0.01). 

Summarizing our results we see that the saturation timescale lies in the range: 
t e {10 8 ; 10 10 }s. A next step is to specify how efficient the twisting of field lines is. 
For this purpose it is relevant to examine an accretion process on AGN, estimate the 
corresponding evolution timescale, and compare with that of the saturation. 

Since the accretion is related to the self gravitating mass, M sg , King & Pringlc 
considering fuelling of AGN, showed in 16J that M sg is given by 



M sg = 2.76 x 10 5 



-2/27 / e x -5/27 / I 



5/27 



(- 



M 



BH 



0.03/ V0.1/ VO.lii 

23/27 

M P) , (27) 

Vl0 8 M Q / 

where rj is the Shakura-Sunyaev viscosity parameter [17j and e is the accretion 
parameter related to the accretion mass rate M and the bolometric luminosity 

6 EE — — . (28) 

Mc 2 

By defining the accretion timescale to be t evo i = M sg /M one can show that [16] 

-2/27 / e \ 22/27 / £ \ ^ 22 / 27 



t evol = 3.53 x 10 13 ( — ) I — ) 
V0.03/ \0.lJ 



0.1L, 



( 



M B H 



vio 8 m ; s - (29) 

From this expression is clear that the evolution timescale varies from ~ 10 12 s, 
(Mg = 1 and L/L E = 1), to - 10 15 s (M 9 = 0.001 and L/L E = 0.01), where 
Mg = M E h/(10 9 x Mo). If one compares these timescales with that of the transition, 
one finds that t evo i exceeds r by many orders of magnitude, illustrating the high 
efficiency of the curvature drift instability. 

We have shown that the saturation of the CDI is extremely efficient, however, it is 
worth noting that the twisting process of magnetic field lines requires a certain amount 
of energy and a natural question arises: is the energy budget enough to provide the 
aforementioned process?. 

Let us introduce the maximum value of the possible AGN luminosity L rnax = Mc 2 
and compare this with the " luminosity" corresponding to the twisting of the magnetic 
field lines L m = AE m /At » AE m /r. AE m is the variation in the magnetic energy 
due to the curvature drift instability. 

We consider AGN with L = 10 45 erg/s, then, assuming e = 0.1, one can show 
from Eq. ([28]) that 

L max = 10 i6 erg/s. (30) 

If the reconstruction of the magnetosphere is feasible, then one has to satisfy the 
condition L m < L max . The magnetic "luminosity" can be estimated straightforwardly 

Lm = §- T Rl*. (31) 
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Figure 5. The dependence of L m / L max versus 8. The set of parameters 
is n = 3 X lu -5 s -1 , L = I0 45 erg/s, 7 e o = {10 5 ; 10 6 ; 10 7 ; 10 8 }, R B « R c , 
n e0 = 0.001cm" 3 , \ e = W0R C and \ x = R c . 

where AV^ ps R^AR = R%k (k = AR/R C « 1) is the volume, where the twisting 
process takes place and B r behaves as 

B r = B° r e t/T , (32) 

where is the initial perturbation of the toroidal component. 

Let us consider the following set of parameters 7 e o = {10 5 ; 10 6 ; 10 7 ; 10 8 }, Rb ~ 
R Cl n eQ = 0.001cm" 3 , A = 100i? c , X x = R c and L = 10 i5 erg/s. In Fig. Owe plot the 
behaviour of L m /L max versus the initial perturbation for the moment of the transition 
(t w t). Different curves correspond to different Lorentz factors. As is clear from Fig. 
[5J the total luminosity budget, exceeds, by many orders of magnitude, the magnetic 
luminosity, indicating that this process is feasible. 

4. Summary 

We summarize the principal steps and conclusions of our study to be: 

(i) We have studied the centrifugally driven curvature drift instability and the 
corresponding saturation process for AGN magnetospheres. 

(ii) By linearizing the Euler, continuity and induction equations we have derived the 
dispersion relation of the parametrically excited curvature drift instability and 
obtained an expression of the growth rate for the light cylinder area. 
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(iii) As a next step we have estimated the transition timescale of quasi-linear 
configuration of magnetic field lines into the Archimedes spiral. Such a shape of 
field lines guarantees the force-free dynamics of particles and necessarily provides 
the saturation of the instability. 

(iv) We have shown that the transition timescale was lower by many orders of 
magnitude than the accretion evolution timescale, illustrating extremely high 
efficiency of the CDI. 

(v) Analyzing the energy budget, it has been shown that the reconstruction of the 
magnetic field lines requires only a tiny fraction of the total energy budget, 
indicating that the CDI is a working process. 

In the present paper we imposed several restrictions. Namely, we described the 
plasma dynamics based on a single particle approach. Although it is natural to suppose 
that the collective phenomena may influence the process of the development of the 
curvature drift instability and a generalization of the present model to a more realistic 
astrophysical scenario is needed. 

In this paper we considered the magnetic field lines located in the equatorial 
plane. However, in realistic astrophysical cases, the magnetic field lines also might 
have poloidal components. Therefore, to extend the approach, it is essential to study 
the CDI for a general configuration of the magnetic field. 
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